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ABSTRACT 



The propagation of simple acoustic transients in an 
isovelocity water layer was investigated for the purpose of 
studying a relationship between the propagation of transients 
in ducts and in layers. The Laplace transform method w as 
used to obtain approximate solutions for the acoustic 
pressure, particle velocity, and particle displacement result- 
ing from a step-function input in velocity and the particle 
velocity resulting from a gated sine-wave input in velocity. 
Computer programs were written to evaluate and graph 
the resulting waveforms. The waveforms resulting from 
step-function and gated sine-wave inputs in velocity were 
observed and compared with the predicted waveforms. With 
the use of Mylar transducers, good correlation was obtained 
in the region of validity of the theoretical solutions. It 
was found that previously published solutions for the 
propagation of transients in ducts became the leading terms 
of the associated solutions for the layer. 



LIBRARY 

NAVAL POSTGRADUATE SCHOOL 
MONTtREY, CALIF. 93940 



TABLE OF CONTENTS 

Section Page 

1. INTRODUCTION 13 

2. THEORY 15 

Development of Solutions 15 

Investigation of Solutions 28 

3. EXPERIMENTAL APPARATUS 41 

Tank 41 

Transducers 1+3 

Electronics Equipment 44 

4 . EXPERIMENTAL INVESTIGATION 50 

Experimental Procedures 50 

Analysis of Data 54 

Error Analysis 6 l 

5. RESULTS AND CONCLUSIONS 71 

6 . BIBLIOGRAPHY 72 



3 



LIST OF TABLES 



Table 

3-1. 

1 + . 1 . 



Page 

Electronics Equipment 47 

Tabulated Errors for Quantities Displayed 70 
in Figures 4*2 - 4*5 



5 



LIST OF ILLUSTRATIONS 



Figure Page 

2.1. An Isovelocity Layer 15 

2.2. Fragment of Radial Particle Velocity 36 
vs. Delayed Time, t‘, as Predicted by 

the Computer Program 

2.3. Fragment of the Radial Particle Displace- 37 

ment, , vs. Delayed Time, t’, as 

Predicted by the Computer Program 

c/a = 0.35 

2.2f« Fragment of Radial Particle Displace- 38 

ment, Jv-n , vs. Delayed Time, t r , as 
Predicted by the Computer Program 
c/aco n = 1.00 

2.5* Fragment of Radial Particle Displace- 39 

ment, ,§ rn > vs. Delayed Time, t*, as 
Predicted by the Computer Program 
c/aco n = 3.00 

3.1. Laboratory Setup 42 

3.2. Geometry of Layer (Top View) 4 2 

3«3* Transmitters and Receivers 45 

3.4. Block Diagram of Electronics Equipment 46 

3-5. Load and Bias Network 49 

4.1. Observed Waveforms 53 

4.2. Experimental and computed £ 7 * vs. 56 

for Particle Displacement Waveform 
Generated by a Step Input in Velocity 

c/aco n = 0.35 

4.3. Experimental and computed VS. 27 ^ 57 

for Particle Displacement Waveform 
Generated by a Step Input in Velocity 

c/ao^, = 1.00 



7 



LIST OF ILLUSTRATIONS (Cont'd) 



Figure 
1*4* * if • 

4-5. 

4.6. 

4.7. 
4*8. 
4.9. 



Experimental and computed vs. 

for Particle Displacement Waveform 
Generated by a Step Input in Velocity 
c/acOr, = 3*00 

Experimental and computed ^/^^v s. c/aw„ 
for Particle Displacement Waveform 
Generated by a Step Input in Velocity 

Ratio of Experimental To Observed 
Normalized Particle Displacement Amplitudes 
vs. Delayed Time Generated by a Step Input 
in Velocity 

Directivity Pattern for 1.0 cm Radius 
Transmitter 

Directivity Pattern for 2.5 cm Radius 
Transmitter 

Equivalent Circuit of Transmitter 



Page 

58 

59 

60 

62 

63 

67 



8 



LIST OF SYMBOLS 



radius of source 



c = c 



free field speed of sound in water 

— J5l 

nth mode 



j group velocity for frequency u> in the 



c = 



n 



c^/c^ , phase velocity for frequency u> in the nth mode 
depth of the isovelocity layer 
mode number 



— P . acoustic pressure 

1 at 5 



r,z,t cylindrical coordinates, assuming radial symmetry, 
and time 

co„ angular cut-off frequency for the nth mode 

^ density of water 

acoustic velocity potential 

P = 

TT = V $ particle velocity 

§> - Z„< z) (r,t), acoustic velocity potential for the nth 

n ° mode 

Z ^ (z), transverse function for the nth mode in the layer 

(r,t), radial portion of velocity potential for the nth 
mode in the layer 

p^ acoustic pressure for the nth mode 

particle velocity for the nth mode 



u 



n 



u 



, radial component of particle velocity 
for the nth mode in the layer 






radial particle displacement for the nth mode 



rn - 



u 



It r\ 



/Z (z)^ fragment of the radial particle- 






velocity for the nth mode in the 
layer 



9 



LIST OF SYMBOLS (Cont'd) 



P 



rn 



fragment of the acoustic pressure for the nth 
mode in the layer 



5 



rn 



fragment of radial particle displacement for the 
nth mode in the layer 



ft n ( z,t) = /R ^ (r) longitudinal portion of the velocity 

potential for the nth mode in the duct 



R n (r) transverse function for the nth mode in the duct 

p = p n / R n (r) fragment of acoustic pressure for 
ir r the nth mode in the duct 



v an = R n (r) fragment of longitudinal particle 

r * velocity for the nth mode in the duct 

fragment of longitudinal particle displace- 






ment for the nth mode in the duct 

longitudinal particle velocity for the nth mode in 
the duct 



K »H 



modified Bessel functions of the first kind 
Bessel functions for the first kind 
modified Bessel functions of the third kind 
U v (w+,ARG) Lommel functions of two variables 

““'Si 

, argument of modified Bessel functions 



.[-? 

ARG= ©s CO n yt -T , argument of Bessel functions 
$( s 



fs 2 + 



T = r-a , time of flight of signal from source at r=a to 
c receiver in layer 

= z/c, time of flight of signal from source at z=0 to 
receiver in duct 



t* = t-T delayed time which represents the elapsed time 
since a signal traveling with velocity c first 
arrived at position r in the layer or z in the duct 

(t) = v rn (a,t), velocity input at r=a in the layer 



10 



LIST OF SYMBOLS (Cont'd) 



(t) = v^ n (0,t), velocity input at z=0 in the duct 

0 (t) = 0 (0,t), pressure boundary condition at z=0 in the 

duct 



1( z) = 0 for z 
= 1 for z 



* °) uni 

> 0 ) 



unit step function 



fW 






o m 



A 



e 

E c 

R 

R, 



s 



e 



relative phase of layer solutions with respect to 
corresponding duct solutions 



rr\ 



values of t' for which S i 
values of t' for which 
area of transducer face 



0 in the layer 
0 in the duct 



blocked capacitance of transducer 

compliance of transducer 

input voltage to transducer 

polarization voltage of transducer 

source resistance in transmitter circuit 

load resistance in the receiver circuit 

equilibrium spacing between transducer film and 
back plate with polarizing voltage applied 

transducer electro-mechanical transfer function 



11 



ACKNOWLEDGEMENTS 



The continuing advice and assistance of Professor 
Alan B. Coppens of the Naval Postgraduate School is 
gratefully acknowledged. 

Thanks are also due to Mr. Ernest G. Janzen for 
his assistance in the preparation of the computer programs. 

This research was supported in part by the Naval 
Ships Systems Command. 



12 



1 . 



INTRODUCTION 



While most investigations of the propagation of energy 
in waveguides have been done for electromagnetic cases, 
some have considered the propagation of pulses of acoustical 
energy in waveguides^ * ^ » 3 » 4) , The dispersion of an ultra- 
sonic pulse propagated through a waveguide of circular cross- 
section filled with water has been used for the approximate 

determination of the frequency spectrum of the original 

( 1) 

pulse . It has been showi that the dispersion of a 
rectangular pulse along an acoustic waveguide of rectangular 
cross-section can be used to calibrate the waveguide so it 
can be used to predict the effects of similar dispersion on 
the transmission of pulses of other shapes that are not 

( p ) 

amenable to actual computation' ' . Pulse-compression 

studies in an acoustic waveguide show that a frequency 

modulated input pulse can be compressed into an output pulse 

of shorter duration and of higher peak-amplitude^). In 

addition to the waveguide experiments, investigations of 

transients in layered media have been made^ . W. C. 

( S ) 

Knudsen' J ' investigated the propagation of a pressure 
transient in a surface layer of liquid which lay over a deep 
bottom layer of liquid. Many of these studies have used 
approximation techniques such as the method of stationary 
phase. A summary of these techniques useful for calculating 
the resultant propagated waveform has been presented by 
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<J . R. Wait(?). The original work of Rayleigh on scatter- 
ing and radiation of harmonic disturbances by small 

cylinders was extended to transient disturbances by <J . W. 

( 8 ) 

Miles using the Laplace transform' R. B. Barakat 

used the Laplace transform to investigate the radiation of 

acoustic pulses of various shapes from an infinitely long 

( Q ) 

circular cylinder v 7 ’ . 

Recently, exact solutions for the propagation of simple 

transients in waveguides (step-function and gated sine-wave 

inputs) have been published^’ ^ ^ . Some work has also 

been done which suggests a correlation between these exact 

solutions in ducts and solutions for the propagation of 

( 12 ) 

acoustic transients in isovelocity layers 

The purpose of this research is to investigate the 
region of validity of the suggested correlation and to see 
what, if any, restrictions must be placed on the resulting 
solutions describing the transient response of the layer. To 
facilitate this comparison, those results of the propagation 
of transients in ducts' * which will be of value are 
summarized in Appendix A. 
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2 . 



THEORY 



Development of Solutions 

The response of an isovelocity layer to simple transient 
excitations will be developed using the following assumptions: 

(a) Cylindrical velocity-source of radius a. 

(b) Constant layer depth, d. 

(c) Perfectly-Reflecting upper and lower boundaries. 







An Isovelocity Layer 
Figure 2 . 1 

For radially-symmetric propagation in the layer indicated 
in Figure 2.1, the wave equation for the acoustic velocity 
potential ^ can be written in cylindrical coordinates as 



( t-' + J- Tr * |r* - ~i §£.) $ = 0 



( 2 . 1 ) 



where r is the radial direction and z is the transverse 
direction. 
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The separation of variables 



£„§„Cv,t} *• Z n i n ()) 0 h ( r> t) ( 2 - 

leads to two partial differencial equations: one for a 
"transverse” function 2^ which depends only on the height 
z of the layer 

J (2- 

and one for a "radial" portion 0 n of the acoustic velocity 
potential, 






(2.4) 



which depends on both time, t, and the distance from the 
origin, r. The value of the separation constant C0 n depends 
on the nature of the boundaries at z = 0 and d. For 
example, if we have a pressure-release top and bottom, 
then p = =0 at z = 0 and d and Eq. 2.3 has the 

solution 






(2.5) 



where O) n = 



n = n7rc /ji * 



Furthermore, since the source is 



assumed symmetric with respect to the plane z = d/2^just 
those modes having pressure antinodes at d/2 will be excited. 
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Thus only odd values of n are allowed. Therefore Eq. 



2.5 can be rewritten as 



( l ™>) ^ n. 1,3^; e,1c.. ( 



2 . 6 ) 



Since this transverse function affects the form of Eq. 2.4 
only through the allowed values of CJ n , it is convenient to 
establish a set of reduced acoustic variables which suppresses 



the function 2 ^ 



The well-known acoustic relations 






C TTJ 



(2.7a) 



^ ^ (2.7b) 

and 

where n is the mode number, p n the acoustic pressure for 
the nth mode, u^ the particle velocity for the nth mode, 
and u rn the radial particle velocity for the nth mode can be 
used with Eq. 2.2 to show that 

■T* = -C ^^t). « 2 - 8 > 
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The "fragment " p^ n of the acoustic pressure for the nth 
mode is defined as 




where 0^ = 0^ (r,t) is the radial portion of the velocity 
potential for the nth mode which was previously defined by 
Eq. 2.2. Combination of Eqs.2.2 and 2.7c allows the 
definition of the "fragment" of the radial particle-velocity 
for the nth mode, 




Since the transverse function (z) is independent of 

time, the radial particle-displacement for the nth mode may 
be obtained by integration of u rn with respect to time. This 
enables the "fragment" of the radial particle-displacement for 
the nth mode to be written as 



J . J 

These reduced acoustic variables or "fragments" 0 n , , 

p , and \ will be used throughout the remainder of this 
rn’ 0 rn 

development . 







X 

LL 



rn 



dt ( 2 . 11 ) 
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Since we are considering the propagation of transients, 
Eq. 2.4» a differential equation in two unknowns, cannot be 
solved by separation of variables. The Laplace transform 
method will be used to obtain the desired solutions. 

The Laplace transform of Eq. 2.4 with quiescent 
initial conditions results in a differential equation in one 
unknown which may be written as 





O 



( a. 12) 



where the bar over the function denotes the Laplace trans- 
form of that function: 



<V r - s > = Xf&T,*)} 1 <PJnt)dt, 

The change of variable 



/*- 



jr 

C 



(2.13) 



casts Eq. 2.12 into the form 



* _ 

>•- 

The solutions to Eq. 2.8 are I 

o 

modified Bessel functions of the 
Conservation of energy requires 




(2.14) 



(/J. ) and K q (/^), the 
first and third kind^^. 

to go to zero as s-^oo , 
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so that the only acceptable solution of this differential 
equation is K q ) . The Laplace transform of 0 n can 

now be expressed as 

% = c'kjsa (2 - 15) 

where C' is an amplitude constant to be determined from 
the boundary conditions. 

With the use of the relationship expressed by Eq. 2. 7c, 

the Laplace transform of AT becomes 

rn 



h r 1/jl hr 




( 2 . 16 ) 



The known boundary condition at the source, located at 
r = a, for propagation in the nth normal mode, 



= V<J*) , 

Sr ’ 

has the transform 



(2.17) 



^rts Vc*. • 



(2.18) 
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Equations 2.16 and 2.18 allow evaluation of C* , 



c - 



-7, 



o_C. 






(2.19) 



so that ir rr can be expressed as 



K,(t\F^) 



( 2 . 20 ) 



The substitution for C 1 in Eq. 2.15 allows to be written 



as 






-t 



k.(iVsw) 



( 2 . 21 ) 






The expression for the Laplace transform of the fragment 
of acoustic pressure, P rn > can be obtained from Eq. 2.21 
if we recall that >f© r ^ ~ ^ ^^2 corresponds to 

in the transform domain. Thus, 



jf rn ^ 

\fsW ; KizU^n') 



( 2 . 22 ) 
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The Laplace transform of the particle-displacement 
fragment, ^ , can be obtained from / V" rn by either 

integration with respect to t in the time domain or multi- 
plication by 1/s in the transform domain. The result is: 



£ 



nr\ 



V. 




s 


K,(tVsW) 



(2.23) 



Inverse transforms of Eq. 2.20 - 2.23 have not been 
tabulated, so approximation methods must be used. 

If ( Eq. 2.13) is greater than unity the asymptotic 

f orm 




(2.24) 



•( ±2 13 ) 

may be used * . The requirement /*> 1 can be 

satisfied by restricting c/au> n to values less than unity. 
Substitution of Eq. 2.21+ into Eq. 2.20 yields 





(2.25) 
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which simplifies to 



^rn'V £Kp£-T|/«V«o] (j + 

\)* 



( 2 . 26 ) 



where T = r-a is the time-of -flight of the signal from 
c 

source to receiver. Since c/a<O n is restricted to values 



less than unity, the denominator can be expressed as 



I 



0 



+ 



& c yfsv.v*' j 



(2.27) 



so that we have 



v r„=Vj(F ^xpL-TVs 1 ^] [i+ i’tfTw’J Li gfy/iwj < 2 - 



28) 



which simplifies to 



C expj^Ti/ — \Za, 3c. (l~y) 



(2.29) 






where terms of order (c/a to ) have been dropped. 

n 
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With similar manipulations the other desired acoustic 



variables can be obtained from Eqs . 2.21 - 2.23: 

( IXpfcWsW] / (+ *\ €X (2.30) 



^ - lc=^-L.(b g) ( 2 . 31 ) 



|F V^cxpL-t/^?] _ 

s 



±L^— l\-2\ w[' t N 



(2.32) 



Equations 2.29 - 2.32 approximate the Laplace transforms 
of the desired acoustic variables for the boundary condition 
V a (t) at r = a with the restriction c/a to 1. 

Comparison of Eq. 2.29 with Eqs. A. 7 and A. 9 shows 
that the terms of Eq. 2.29 have known inversions. For a 

i 

step input, for which V & = /g the inverse transform for 
the fragment of the radial particle-velocity can be written 
down directly from the known results of the duct (see 
Appendix A. 
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■where ARG ~ ^ n \j t^-T^ and it is always understood that 
all solutions vanish for t £ T. 

The inverse transforms of the other acoustic quantities 
can be found similarly by using the known transform pairs 
from the duct case and the additional transform pairs 
developed in Appendix A. The results are 




_ c. 






u 2,(z*-t) 



(2.34) 




(2.35) 




(2.36) 



The quantity t* = t - T appearing in Eqs . 2.33 - 2.36 and 

2 2 

in ARG as ( t - T ) ( t 4 T ) = t - T represents the time 
which has elapsed since the signal traveling with velocity c 
first arrived at the position r. It will be referred to as 
the delayed time. 

Equations 2.33 - 2.36 are the solutions for the transient 
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response of an isovelocity layer to a step input in velocity 



at r = a and t = 0 . 

Similar procedures may be used to obtain the fragment 

of the radial particle-velocity for a gated sine- wave input 

77 0-5 

in velocity for which Y^- * 




( 2 . 37 ) 



where = Lommel function of two arguments^ , 



and 




'vxr±® 



Other acoustic quantities could not be completely solved in 
this case. Details are shown in Section A-3 of Appendix A. 

Comparison of the above solutions with the solutions 
for the duct listed in Appendix A shows that in all cases 
the solutions previously obtained for the duct become the 
leading terms in the solutions developed for the layer. 

A physical understanding of this can be gained by 
comparing propagation processes in a duct with those in a 
layer. While a duct presents a constant cross-section to 
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the propagating wave form for all time, a transient 
propagating outward in a layer is presented with an 
effective cross-sectional area which increases as r. How- 
ever, for sufficiently large r the cross-sectional area will 
not change appreciably as the waveform propagates through 
unit distance. In the limit, therefore, we expect the 
radial divergence contained by the term” in Eq. 2.4 

to become insignificant. In this limit, propagation processes 
can be described by an equation virtually the same as Eq. A3 
which describes propagation in a duct. 

For this reason it is not surprising that the solutions 
for the duct become the leading terms of the layer solutions. 
The remaining terms in the layer solutions are correction 
terms necessitated by the fact that the layer does not 
actually present a constant cross-section to the propagating 
wave form. 

The solutions for the transient response of an 
isovelocity layer have been developed with the restriction 
c/a (2> <. 1. Since the cutoff frequency is given by 

n Tr C./j , this restriction places a limit on the 
relation between the depth of the layer and the size of the 
source, d/a y (oTT • If it is desired to remove this restric- 
tion a different expansion of the ^ <yi 's involved must be 
used. 
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Investigation of Solutions 



The solutions developed for the layer have been obtained 

f or a single mode. In determining the contributions to the 

solutions from other modes, attention must be given to 

the transverse behavior of the transmitted and received 

signals. The boundary condition on velocity at the source 

at r = a says that V (t) vanishes for z <.0 and z >d and 

is constant for 0 < z < d. This condition can be met by 

representing V (t) by a spatial Fourier-superposition of 
a 

jjl(z) - U(z-d)] in terms of the normal modes across the 
face of the source: 



V» ■ 4 Ui-) 



S/n/2r>-/)7T%. 

% 



( 2 . 38 ) 



Notice that the higher modes contribute as for this 

input. These higher modes have positive and negative 
portions in their waveforms over the depth of the layer. 
These positive and negative portions will partially cancel 

it -m k 

over the face of a receiver which extends over the layer 
depth. In effect the receiver will average the received 
signal over the area of its face^"^: 






( 2 . 39 ) 
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Evaluation of the above integral gives 







(2.40) 



Therefore each mode contributes as (^Ti) or 1> 1/9 , 

1/25, l/49, for the first four excited modes. 

A computer program was written to evaluate the 
effects of these higher modes on the solutions developed. 
Computer predictions considering (a) the fundamental mode 
alone and (b) the fundamental mode plus the next four non- 
zero higher modes, were plotted versus delayed-time . 
Comparison of these two predictions as presented graphically 
by the computer showed changes in the amplitude of less 
than 1 per cent and shifts in axis crossings of less than 
0.01 local wavelengths. Since these effects are negligible 
compared to the experimental errors encountered in this 
research, it was sufficient to retain only the fundamental 
mode ( n = 1 ) . 

Examination of Eq. 2.36 shows that it consists of 
two terms. The first (main) term is from the duct 
solution and the second (correction) term is from the layer 
effect. An estimate of the effect of the correction term 
can be made by combining the leading terms of the indicated 
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sums. Equation 2.36 can be rewritten as 




(2.41) 



Under the assumption t' = t-T « T, which limits considera- 
tion to the leading portion of the transient waveform, we 
have t^T, and Eq. 2.41 becomes 



^ Srn^ 1 3*1 W (2.42) 

where © = ARG = t' 2.T 

and ett) = ( '■ %) fc) 4 • 

Substitution of trigonometric approximations of the Bessel 
( 13 ) 

functions in Eq. 2.42 gives 



jr 




* 




cos(©-lIT') - c.cs($-£?y (2.43) 



By combining the cosine terms of Eq. 2.43> we find the 
expression for the fragment of radial particle displacement to 
be 



\J 



ir 1 

7* «- 



to 



*T 




costin' + ^ 



(2.44) 
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where 



= tan^B ( t ) . 

Equation 2.1+1+ shows that the correction term in the 
layer solution introduces a phase angle that was not 

present in the duct solution. The zeroes of the right hand 
side of Eq, 2. 1+4 give the times for which the particle 
displacement vanishes. It has been shown that the duct 
solution becomes the leading term of the layer solution; 
therefore, the zeroes of the first term of Eq. 2.1+3 give the 
nulls, or axis crossings, of the particle displacement in a 
duct. If these times are evaluated, the shift in time of 
associated axis crossings between the layer and the duct 
solutions can be predicted. 

The zeroes of the layer solution, designated by » 

are obtained when 

COS (6-2E+ ^ = 0 > (2.1+5) 

or 

- (^^') Tr 3 m = l ^^ ctc * (2 * 46) 
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The zeroes of the duct solution 'will be designated by ^ 
They are obtained when 



COS 





( 2 . 47 ) 



or when 




2, e'te. 



( 2 . 48 ) 



Since we have t' < ^ T, c/a co <. 1, and r ZZr a, so that 



-w x - m = <x i ^ , 2 .49) 

a small-angle approximation can be applied with the result 

+o- n * = s' = r (2.5o) 



where 




2.-% 

h ^ r > yfpr 7 
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into Eq. 2.2+6 and 



Substituation of this expression for 
evaluation of © in Eqs. 2.2+6 and 2.2+8 gives 



(co^/srF -h H ) 




( 



Zm- I 

Z 




(2.51) 



and 



-( z ¥^) 7r 



(2.52) 



Manipulation of Eqs. 2.51 and 2.52 reveals that the ratio 



predicting the shif t in axis crossings. 

This result is approximate, since only the first terms 
of the sums in Eq. 2.36 were used to develop it. To 
evaluate the accuracy of this approximate result, computer 
programs were designed to evaluate the predicted waveforms 
of Eqs. 2.33 and 2.35 - 2.37* They were written so that 
the main terms and the correction terms of these equations 
could be evaluated separately. This permitted comparison of 




can be expressed as 







(2.53) 




is a convenient parameter to use in 
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the layer solutions and the duct solutions and showed the 



behavior of the correction term as a function of c/a W 

n 



These programs, one for the step input (Appendix A) and 

one for the gated sine-wave input (Appendix D) consider 

up to 11 terms of the sums involved. The programs give 

printed and graphical output of the main terms, the 

correction terms, and the complete expressions of the 

equations mentioned above. Figure 2.2 shows the velocity 

fragment predicted by the computer. Figures 2.3 - 2.5 

show the complete expression, the main term, and the 

correction term of the displacement fragment for several 

values of c/aiJ 

' n 

Information about the envelope of the displacement 
fragment waveform can also be obtained from Eq. 2.44 
which gives the amplitude as \| l "V 6 Ct) . Under the 
assumption t' << T, B(t) becomes 




(2.54) 



and 




(2.55) 
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Therefore, the amplitude of the initial portion of the 
waveform under consideration is 




The approximate solutions for the propagation of simple 
transients in an isovelocity layer have been obtained, and 
computer programs to evaluate these solutions and predict 
the axis-crossings and amplitude-modulation in the layer 
have been written. 
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(b) Gated sin* wav* input with c/*h,= 1»00 

cm, a • 2»5^ cm, d*7»97 cm* <•> • kHz< 

Fragment of radial particl* velocity, v^ , va. delayed time, t 1 . 



Figure 2.2 
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(b) Layer correction term 

Fragment of the radial particle displacement, § rn , vs. delayed time, 
t', as predicted by the computer program, c/aco^r r-6l.5 cm, 

a = 2.54 cm, d- 2.80 cm. 



Figure 2.5 
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Fragment or the radial particle displacement, 5«-n, vc» delayed time, 
t', ac predicted by the computer program. c/ao H =1.00, r- 61.5 cm, 
a •=2.5^ cm, d-7*97 cm. 



Figure 2.4 
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Conpleco or.pre oion (layer solution) 
Main tarn (duct aolution) - ~ ' 



Fragment of the radial particlo displacement, Sr*> vs * delayed time, 
t‘, as predicted by the computer program. c/au>,= 3*00j r=6l .5 cm, 
a- 1.00 cm, d = 9*92 cm. 



Figure 2.5& 
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( Arbitrary Units) 




Layer correction tern 



Fragment of the radial particle diaplacement,^ , vs. delayed time, 
t', as predicted by the computer program. c/aw s =5.00, r - 6 I .5 cm, 
a = 1.00 cm, d = 9«^2 cm. 



Figure 2 . 5 b 



2+0 



3. 



EXPERIMENTAL APPARATUS 

The experimental apparatus used in this research 
consisted of a tank, several transmitters and receivers, 
and the associated electronics equipment. The laboratory 
setup is shown in Fig. 3 . 1 . 

T ank 

The water layer was contained in a tank constructed 
of 3/4 inch marine plywood which was 30 cm high, 1 m wide, 
and 1.5 rn long. The plywood was coated with Varathene 
to provide waterproofing, and lined with 1/4 inch 
polyethylene foam. This provided the desired pressure- 
release bottom boundary for the layer. The water-air 
interface provided the top pressure release surface. With 
the sides and ends of the tank also lined with polyethylene, 
the tank could be excited as a pressure-release cavity which 
served to accurately determine the depth of the layer. 

This technique is described in more detail in section 4* 

Cylindrical sources could not be used because the tank 
was too small. Transmitters were therefore built in the 
shape of semi-cylinders and mounted on a rigid plate. A 
neoprene-coated steel plate, 3/8 inch thick, was made to 
fit in one end of the tank to act as the rigid boundary. 

This rigid surface served as an acoustical mirror, so that 
the reflected pressure appeared to be emitted from an 
image source (an identical semi-cylinder back-to-back with 
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the real source) . Under this condition cylindrical symmetry 
was approximated. Figure 3*2 illustrates the geometry 
of the layer. 

Transducers 

Mylar dielectric transducers were selected to be used 
in the layer for the following reasons? 

(a) Their broad-band transmitting and receiving 
characteristics allowed the transmission and reception of 
the transient excitations without introducing significent 
frequency-sensitive modulation of the signals. 

(b) They could be conveniently constructed to the 
required semi-cylindrical shapes and sizes. 

(c) Previous work^ has indicated that, when 

driven with a low-impedance voltage source, they generate 
an acoustic signal for which the normal particle velocity at 
the source is proportional to the input voltage. 

Two transmitters were constructed to the required 

( 17 ) 

semi-cylindrical shape, using the techniques described in , 
one with radius 1.0 cm and one with radius 2.5 cm. Both 
transmitters were 12 cm high which corresponded to the 
maximum layer depth to be encountered. These dimensions 
were chosen in order to cover the desired range of values 
for c/a W n . 

A Mylar dielectric receiver, identical in nature to the 
transmitters, but with a plane, rectangular face measuring 
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3 cm x 1 4 cm, was used to measure the averaged (over 
depth) radical particle velocity and displacement in the 
layer. By proper adjustment of its electrical load, the 
receiver could be made sensitive to either particle velocity 
or displacement. When terminated by the 1 megohm 
impedance of the voltage amplifier, it was displacement 
sensitive; when shunted with an impedance of 100 ohms, it 
became sensitive to velocity. Transmitters and receivers 
are shown in Figure 3*3* 

Two pressure-sensitive receivers were constructed, 
but neither proved to be satisfactory. A probe using a 
1/8 inch barium titanate cylinder as the sensing element 
was built, but its low sensitivity resulted in an unacceptable 
signal to noise ratio. A 2.5 cm x 2.5 cm sensing element 
cut from a 0.5 cm Glennite disc was also constructed. 
However, the resonant frequency of the element was 
sufficiently close to the frequencies at the front of the 
received waveform so that ringing of the element caused 
severe distortion of the received signal. For this reason, 
and because it was felt that the velocity and displacement 
measurements would adequately describe the waveforms, 
pressure measurements were not taken. 

Electronics Equipment 

Figure 3*4 shows the block diagram of the equipment; 
Table 3.1 identifies the items in the figure. The oscillator 
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un the left, a t.p cn radius baclrplate and. finished 
transmitter; center, partially completed and finished 
receivers; on the ri~ht, the l.J era radius transmitter. 

transmitters and ccoivers 

figure d • p 
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Figure 3.4 



Table 3.1 





Electronics Equipment 


osc 


Audio oscillator, TS-382 D/U 


TBG 


Tone burst generator, General Radio type 1396-A 


UPG 


Unit pulse generator, General Radio type 1217-A 


FREQ 


Electronic counter, Hewlett Packard 5233-!-' 



PWR AMPPower amplifier, Hewlett Packard 46 7 -A 



AMP 1 


Voltage amplifier, Teletronix type 1121 


AMP .2 


Wide band chain amplifier, SKL model 202-D 


RC FIL 


Homemade variable low-pass RC filter with 
400 kH z roll-over. 



SKL FIL Electronic band pass filter SKL model 302 



SCOPE 


Dual beam oscilloscope, Fairchild model 766 -H 


LBN 


Load and bias network (made for earlier work) 
shown in Fig. 3*5 
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and tone burst generator were used as the initial signal 
source for the gated sine-wave input. For the step input, 
the unit pulse generator was used. 

The received signal was sent from the receiver 
through the load and bias network, amplifiers, and filters 
to trace "B" of the dual beam oscilloscope. The input 
signal taken from the output of the power amplifier was 
displayed on trace "A". Measurements of the received 
waveform axis crossings and amplitudes were made directly 
from the oscilloscope display, or photographed for later 
analysis. Triggering of the scope was provided by gating 
signals from either the tone burst generator or the unit 
pulse generator. The load and bias network shown in 
Figure 3.5 served to supply the. 300 volt bias needed by 
the Mylar transducers without allowing this voltage to 
appear elsewhere in the circuit, and to provide the proper 
electrical load for the receiver so that either particle 
displacement or velocity could be observed. The 22 ohm 
resistor in the transmitter circuit was needed to suppress 
parasitic oscillations caused by the transient electrical 
excitation of the capacitance of the transmitter and the 
residual inductance of the power amplifier. 

The rise time of the system was experimentally 
determined to be about 0.5 M sec. 
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Figure 3.5 



4. EXPERIMENTAL INVESTIGATION 

The purpose of the experimental investigation was to 
(a) observe and analyze the received radial particle-velocity 
and displacement waveforms in an isovelocity layer; (b) 
compare these waveforms with theoretical waveforms produced 
as output of computer programs based on the theoretical 
developments; and (c) vary the parameter c/a £j over a 
range of values to see if the solutions broke down when 
the basic restriction c/a <,1 was violated. The greater 
part of the experimental investigation was made using the 
step-function input; the gated-sine wave input was used 
only infrequently. Experiments were conducted using the 
procedure described in the following paragraphs. Five basic 
variables characterized each experiment; the input function 
(gated sind-wave or step-function), the range r, the 
transmitter radius a, the layer depth d, and the speed of 
sound c. 

Experimental Procedure 

Since c/a W was a fundamental parameter of the 
solutions, the input function and the value of c/a w n to 
be used were first chosen for each run. The 2.5 cm 
transmitter was used for values of c/a 60 n from 0.35 to 
1.00; the 1.0 cm transmitter was used for values ranging 
from 1.00 to 3*00. The temperature of the water was 

( 18 ) 

measured to determine the speed of sound c in the layer' '. 
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With a, c, and the desired value of c/a CO determined 
the necessary value of could be calculated. The tank 

was then filled to a convenient depth and excited as a 
pressure-release cavity of known length and width. The 

depth of water was adjusted until the maximum resonance 

\ 

peak of the fundamental mode was observed at the 
desired . From and the known length and width of 
the tank the effective depth was calculated for the computer 
program which was written with depth rather than (jb n as 
an input. The end plate and transmitter were placed in 
the tank, an input signal was fed to the transmitter, and 
the receiver was placed at the range of 61.5 cm (used in 
all experiments) by measuring the time-of-f light of the 
signal from transmitter to receiver. 

The received waveforms were photographed from the 
oscilloscope display. For the step input, several photo- 
graphs were taken of each waveform. An expanded time 
scale was used for photographs of the leading portion of 
waveform where the high frequencies were contained; a 
more compressed time scale was used for suceeding photo- 
graphs. For a typical experiment, three photographs of 
each waveform were made utilizing, for example, 5 /W.sec/ 
cm, 10 /u. sec/cm, and 50 MSec/cm time scales. For the 
step input experiments, both received particle displacement 
and velocity were photographed in most cases. The photo- 
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graphs were then read using dividers and a scale to determine 
the values of delayed time at which axis crossings 

occured. The amplitudes of the peaks in the waveform were 

■- ! I 

measured in a similiar fashion. Sample photographs are 
shown in Fig. l+.l. Data read from all photographs are 
contained in Appendix B. 

Fbr each experimental run, a corresponding computer 

run was made using the same values for r, a, c, and d 

' * /* 

(and hence c/a W ) as inputs to the computer program. 

As has been shown in Section 2, a solution for particle 
displacement or velocity in the layer can be expressed as 
the sum of the corresponding duct solution and a "layer 
correction." The computer programs could therefore be 
written so that their output for each run contained a 
series of graphs of the duct solution, layer solution, and 
"layer correction" waveforms. The program for the step 
input produced a series of graphs for both particle dis- 
placement and velocity waveforms. The program for the 
gated- sine inpu 1 t'/graphed:"the velaci-tyhseries only. Comparison 
of graphs within a series provided convenient visual evidence 
of the relative importance of the "layer correction" with 
respect to the duct solutions as a function of the parameter 
c/a LO n * A portion of the step input program was used 
to find the axis crossings of the theoretical displacement 
waveforms for the duct and layer solutions. 



52 





steo Input 
in /elocity 

c/au A =o. 

r = 31. cm 

a = . r ; c .'i 

I - 7 . 3 cm 
teirro = * 




late I ine avc 
Input in /elocity 
w= 5 rj . j '-z r Z 

c/ = ] . 

r - - 1 • cn 

o - . cm 

t .,mr> = 6 . 






Another portion of the program evaluated the magnitude 

\ 1 

and delayed time of the amplitude peaks in the layer-solution 

displacement waveform. In this manner the theoretical values 
of ^ and > and the magnitudes of times of the peak 

amplitudes could then be compared with the observed values 
of l m and amplitude peaks found in the experiment with 
the corresponding value of c/a CO . 

Analysis of Data 

Analysis of the response of the layer to a step- 
function input proved to be more involved than was originally 
estimated. For this reason, and because of the limited 
time available for this investigation, data from the runs 
performed with a gated sine-wave input were not analyzed. 
The following analysis pertains to the step-function input 
data only. 

Recall that the simple theory in Section 2 based on 
the trigometric approximations for the Bessel functions 
suggested that if values of 2"^ are plotted versus '%'orr\ > a 
straight line should result. This relationship would offer 
advantages in data reduction because a family of axis 
crossings could be summarized by the slope of a straight 
line. For this reason, graphs which plotted the observed 
and computed values of 2^ versus the computed values of 
were made for each run. Figures 1+.2 - k-k are 
representative examples of these graphs. The solid line of 
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unit slope represents the computed values of • The 

circles are the computed values of , an d the crosses are 

the observed values of . It is clear that each set of 

points can be joined by a straight line which quantitatively 
describes by its slope the predicted and observed shift in time 
of the layer solution axis crossings as compared with those of 
the duct. 

A graph was then made which plotted the predicted and 

^ , 

observed values of * versus c/a Co n (Fig. 4*5)* The simple 

< 2 * 2 

theory would suggest plotting —22 versus a/r (c/aqj ) > how- 






n 



ever, r was fixed for all experiments and only two values of 
a were used. Furthermore it is the validity of the layer 
solutions for various values of c/aco n which is the topic of in- 
vestigation here, rather than the validity of Eq. 2.53* 

Analysis of the amplitude peak data was also made for 
step input solutions. The experimental amplitude peaks read 
from each waveform photograph were first normalized to the 
value of the highest peak in the waveform. The predicted 
peaks were normalized in the course of the computer program. 
The ratio of these normalized quantities was then plotted as 
a function of the delayed time of the predicted amplitude 
peaks for each run. Exact agreement in the shape of the 
transient waveform between experiment and theory would 
result in a line of zero slope. Figure 4*6 shows plots for 
three values of c/att*n” 
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c/a u)k = 0.35 a = 2.5 cm 

r = 61.5 cm d = 2.80 cm 

temp = 20* C • 

Experimental and Computed 7** vs. 7 m* Particle Displace- 
ment Waveform Generated by a Step Input in Velocity. 

Figure 4.2 
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c/atJn = 1.00 a = 1.0 cm 

r = 61.5 cm d = 3.14 cm 

temp = 20 C. 

Experimental and Computed "K* Particle Displace 

ment vv'aveform Generated by a Step Input in Velocity. 

Figure 4.5 
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c / a Wii — 3*00 a * 1.0 cm 

r = 61.5 cm d = 9»42 cm 

temp = 20* G. 

Experimental and Computed 7^, vs. 71m Particle Displace 
ment Waveform Generated "by a Step Input in Velocity. 

Figure 4.4 
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Ratio of Experimental to Observed Normalized Particle Displacement 
Amplitudes vs. Delayed Time Generated by a Step Input in Velocity 

Figure 4.6 



Error Analysis 



> * 

An analysis of the assumption that the semi-cylindrical 
source and steel backplate combination could be used in lieu 
of a cylindrical source was necessary. The complex 
impedance of the plate was evaluated from relations in for 
the range of frequencies contained in the received wave- 

I 

forms. The impedance mis-match was sufficiently large to 
justify assuming that the plate was rigid. The directivity 
of each transmitter when mounted on the steel plate was 
investigated in order to compare the resulting radiation 
pattern with the uniform radiation assumed in the theoretical 
deveopment. Directivity patterns are shown in Figs. 4*7 
and 4 . 8 . In order to determine the effect of the deviation 
of the actual patterns from the theoretical patterns ^ each 
transmitter was clamped with foam rubber so that only a 
90 degree arc in the center of the transmitter was radiating 
sound. Except for small change s in amplitude, no 
differences were observed in the received waveforms as 
compared to those received for radiation over the full 180 
degrees of the transmitter. It was therefore assumed 
that radiation from the transmitters could be assumed 
uniform and that the transmitter-plate combination could 
be used to approximate a cylindrical source. 
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Directivity Pattern for 1.0 cm radius Transmitter 

Figure 4.7 




63 



Directivity Pattern for 2.5 cm radius Transmitter 

Figure 4.8 



The rise rime cf the experimental system was determined 



to be about 0.5 M s © c * Therefore, instead of being a perfect 
step-function input, the velocity at the source was initially 
described by 

_ ex ± 

Vjt) - i - * 

6 —I 

where - 2 x 10° sec “ is the reciprocal of the rise 

time. 

The group velocity of the waveform is given by 

Under the condition <*>>, , and with T) Eq. 4.1 

can be written as 



~r 

7VH* 




By rearranging Eq. 4«J-» the value of delayed time at which 
a frequency co 'will appear in the received waveform can be 
expressed as 



t 




j* 

2. 




t 



It would be expected that the amplitudes of frequencies 
greater than ©<4 conrained in the initial portion of the 
received waveforms would be depleted up to a time given 




and that amplitude distortion of the waveforms would be 
significant up to this time. However, substitution of the 
observed o<f and the largest value of to n experimentally 
encountered into Eq„ 4*2, shows that distortion would be 
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significant at the most for only 1.5 /■ JLsec . Therefore the 
rise time of the system should have no great effect on 
the received waveforms. 

It was necessary to investigate the characteristics of 

the transmitters to see if the velocity of the face 

remained proportional to the step input voltage over the 

interval of time that the received waveforms were studied. 

For times much greater than the rise time, the mass of the 

diaphram can be ignored provided its velocity remains 

reasonably constant over the time interval of interest. The 

resulting equivalent circuit of the transmitter with a source 

( 17 ) 

resistance added is shown in Fig. !+. 9« Writing loop 

equations for the equivalent circuit gives 

e - i- R s - Z +£ jf Jt = » 

which enables us to write the Laplace transform of the 
diaphram velocity as 



^ e. 



Vl- 



(4.3) 



is determined from the electro-mechanical relationships 

(17) 



of the transmitter and is defined as 

c* e 6 
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Nominal equivalent circuit values for the 1 cm radius 
transmitter with a 20 volt seep input applied are shown 
in Fig. 4* 9- Substituting these values into Eq, 4.3 
and taking the inverse transf orm yields 

\f _ o ,1 -** -3^0*. 

- 3.^1 x e 

Since the mass of the transmitter diaphram is ignored in 
this approximation, the particle velocity is a maximum 
immediately after the application of the step input and then 
decays exponentially. Substituting t“360 yuisec (the longest in- 
terval of delayed time over which received waveforms were 
observed) into Eq. 4.3 shows that the velocity has 
decayed to 88 percent of its maximum value at the end 
of this time. The tailing off of the amplitude peaks of 
the received waveforms, as shown in Fig. l±.6, is consistent 
with the decay in velocity of the transmitter diaphram. 

Uncertainties in a, c, r , , and iT^ were the 

possible sources of error in this investigation. The effective 
radius of each transmitter was measured to the nearest 
1/2 mm, which could have introduced an uncertainty of + 2 
percent in the value of a for the 2.5 cm radius transmitter, 
and + 5 percent in the value of a for the 1.0 cm radius 
transmitter. Measurement: of the temperature of water 
to the nearest degree restricted the (systematic) error in 
c to within + 0.5 percent. The method of determining T, 
the time-of -flight of the signal from transmitter to 
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R S = 22 ohm 
C 0 = 1000 pf 
A = 3.8 x 10"3 
.Ac = 1.5 x 10 6 kg/sec-m2 

c s = io ~ 6 m/N 

E Q = 300 volts 

x 0 = 3 x 10"5 m 

p = 10"2 N/volt 

e(.tj = 20 volt step input 

Equivalent circuit of transmitter 
Figure 4.9 
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receiver, gave a possible systematic error of + 0.5 percent, 
which in turn resulted in an estimated error in r of + 0.5 
percent. 

Two sources of error were possible in the determina- 
tion offc) n by excitation of the tank as a pressure-release 
cavity. Observed irregularities of + 1/4 mm in the bottom 
surface of the tank could introduce an error in the local 
effective depth of the layer. Since experiments were 
always conducted in the sarnie portion of the tank, a 
systematic error in UJ would result. This error could 
range from + 0.3 percent' (for the 9*42 cm depth) to + 1 
percent (for the 2.80 cm depth). The random error 
possible in determining the frequency of the resonance peaks 
of the standing wave spectrum proved to be insignificant 
relative to the possible systematic error. 

The experimental values of could be influenced by 
systematic errors on the order of + 1 percent caused by 

^ *1 ), 

nonlinearity and miscalibration of the oscilloscope time-base. 
As mentioned previously in this section, axis crossing times 
were read from a series of photographs with different 
time scales. This method was used so that the initial 
portion of the received signal would always be at the left- 
hand edge of the picture and the axis crossings of interest 
would generally occur in the right half plane of the display. 
This procedure resulted in a random error of + 0.5 percent 
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on all time scales since the display could be read to within 
1/4 scale division out of 50 divisions. 

The errors in the computed values of 2"”^ were influenced 
predominately by the error in to n , since CO n appears 
in the argument of the Bessel functions ( ) from 

which the *27^ ' s are com P u ted. The values of were 

determined when the argument equaled any one of a set of 
constants which resulted in a zero value for the 

amplitude of the particle-displacement solution. Thus the 
argument can be written as 

and if we have tv, «.t , then 

\forr "Cm z ^ ^ 

and 

Therefore the + 1 percent error in (0 ^ could lead to + 2 
percent error in the computed values of r m • 

The errors in the slopes of the curves fitted to the 

> 

experimental 2’U 1 ,S * n Figs ♦ 4*2 - 4*4 were determined mainly 
by the systematic errors in these values of because 

the random errors were minimized by fitting a straight line 

J V 

through the points. Therefore the error in each slope 
was estimated to be + 1 percent. There was a systematic 
error of + 2 percent in the computed values of ; hence 

the error in the slope of the computed 2^ versus *2^^, 
lines was + 2 percent. 
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The errors in c, a, and uJ being independent can be 
combined by the root-mean-square law to yield the error 
in c/a bO n< The percent error in c/ato n varied according 
to the values of a and involved. Table 4«1 shows the 
result of this combination for representative values of 
c/a . Error flags on Figs. 1+.2 -4*5 indicate the 
relative magnitudes of errors in the experimental and 
computed values of and slopes, and the relative 

magnitude of error in c/a oJ . 



T able 4 . 1 



Tabulated Errors for Quantities Displayed on Figures 4 .2-4 *5 





Percent Error 


c/a<o n 


c 


a 


tO n 


c/aco n 


Exp . 
slope 


Comp . 
slope 


0.35 


+0.5 


+2.0 


+ 1.0 


+2.3 


+ 1.0 


+2.0 


1.00 


0.5 


2.0 


0.3 


2.1 


1.0 


2.0 


1.00 


0.5 


5.0 


0.8 


5.1 


1.0 


2.0 


3.00 


0.5 


5-0 


0.3 


5.0 


1.0 


2.0 
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5 . 



RESULTS AND CONCLUSIONS 



The investigation has shown that for c/a 0G n <^l, 
the predicted and experimentally observed shift in axis 
crossings agreed to within + 2 percent. Where the 
restriction c/au> n < 1 was violated, the theory failed to 
predict the observed shift in axis crossings. As c/a ui n 
was increased the failure of the theory became more 
pronounced. 

The correlation between the predicted and observed 
amplitudes was less favorable. The divergence between 
the predicted and observed amplitudes in the initial portion 
of the received signal was unexplained. The decrease in 
the observed amplitudes in the latter portion of each 
received signal is not inconsistent with the roll-off of the 
step-function input in velocity. 

Therefore, with the exception of amplitudes, we can 
say that the theoretical solutions correctly described the 
propagation of simple transients in an isovelocity layer with 
perfectly-reflecting boundaries as long as the parameter 
c/a to n was — 1; and that, although the theory failed for 
values of c/aui n > 1, it did correctly predict the trend 
that the observed waveforms followed. 
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APPENDIX A 



A-l : Summary of the Propagation of Transients in Waveguides . 

Results for the propagation of simple transients ( step- 
function and gated sine-wave input) in ducts which have 
appeared in the literature ^’^) and which will be of value 
in this research are presented below along with some addi- 
tional Laplace transform pairs that are needed. 




Figure A. 1 

Assume a waveguide of uniform cross-section, infinite 
length, and perfectly-reflecting boundaries, which is filled 
with a lossless fluid and excited by a velocity souce at z = 0 . 
See Fig. A.l The wave equation for the acoustic velocity 
potential (D can be written in cylindrical coordinates as 









(A.l) 
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where r is the radial direction, z is the longitudinal direction, 
and radially symmetric processes have been assumed. 

The separation of variables 

= Yl n L(-) <A.2> 

leads to the solution R = J ( 60 r/c) where <J represents 

n o n 

the zeroeth order Bessel function of the first kind, and 
the partial differential equation 

[% -(tT- - G (A.3» 

which describes the propagation in the z direction. 

The value of the separation constant depends on the 

nature of the boundary at r = a (the wall of the duct) and 
the cross-sectional area of the tube. Because the transverse 
properties of the duct affect only 6 j-> and not the form of 
the solution of Eq. A.. 3, all functions of r have been 
suppressed. They can be reintroduced with no difficulty 
when desired. 

The Laplace transform of Eq. A. 3 with quiescent initial 
conditions is 

1 0 (A. k) 

V 
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where the bar over the function indicates the Laplace 



transform of that function, and 0, is a fragment of the 
velocity potential containing z and t. Solving Eq. A. 4 with 
the boundary condition (0,t) gives 

<f>^ ~ (j. ^ exP (jT (A * 5) 

where T = z/c is the time of flight of the signal from the 
source to the receiver. 

The relation CC- n - Rn is used to obtain the 

* *F 

Laplace transform of the fragment of the longitudinal 
particle velocity: 

^kr, =- (A.6) 

c 

Application of the boundary condition on velocity at the 
source v n (0,t) = V o (t) casts Eq. A. 6 into the following form 

- V 0 (A. 7 ) 

The relation 




(A. 8) 
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is evident from examination of Eqs. A. 6 and A. 7. If 
this value of 0 Q is substituted into Eq. A. 5 , the Laplace 
transform of 0 n in terms of the velocity input becomes 

^ = ~clC £*p[-TVsW] (a .c 

The fragment of the longitudinal acoustic pressure is 
obtained from the velocity potential by means of the 
expression n = •'fc j ^ ^ = — ^ • This corresponds to 

P zn = - ^ s0 n in the Laplace transform domain. We 

therefore have the result 




exp 






(A. 10) 



The fragment of longitudinal particle displacement may 
be obtained either by integrating the fragment of longitudinal 
particle velocity in the time domain or by multiplying by 1/s 
in the Laplace transform domain. With the latter alternative, 
the particle displacement fragment can be written as 



§ - \/ e _ (A. 11) 
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Equations A. 7 and A. 9 - A. 11 are the Laplace transforms of 
the acoustic variables of interest. The inverse transforms 



for a step-function input (V Q = l/s) and a gated sine 
— CL? 

input (V = s up- ) are as follows: 

For a step-function input 




c-O 



a * (tzf? 



T c 

— '2»v-| 










Where ARG 




For a gated sine-wave input: 



(rux-j A-(?6)J 

U,(n*r tJ ftefc') + 



wave 



(A. 12) 



(A. 13) 



(A. 14 ) 



(A. 15) 



(A. 16) 



(A. 17) 
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(A. 18) 





(A. 19) 



Where ARG = U> n \J t 2 -T 2 

aov = ui U-t)(/ ^ C a/c > ) 

The U's represent Lomrael functions of two 
variables 




A-2 : Development of Additional Laplace Transform Pairs 

for Step Input in Velocity in an Isovelocity Layer . 

The transform pairs developed in this section are 
those needed for the layer solution in addition to those 
summarized in section A-l. F denotes the transform under 
consideration; f(t) is the inverse transform of F. All f(t)'s 
are defined as zero for t ^ T. 

1 . 

? - expC-Tt/^W^) (A. 20) 
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With V o = l/s, Eq. A. 20 becomes 



F - exp 



(A. 21) 



Equation A. 21 is of the same form as Eq. A. 7 with 

— tJ 



V o 



and tO r . Therefore the inverse 



transform of Eq. A. 7 may be used with C«0 = 
inverse transform of Eq. A. 7 is 

f 0^^-, Ae*) 



The 



(A. 22) 



where ARG = 60> n \f t 2 -T 2 



and 'UJV =. CO ( -£ -T') ^ / ± C j/c^J 



The quantity c g“/c can be written as \[ | - 
With tO — , we have 'Wy. =- njX, - cO n ( t-T ) . 

This in turn means that both terms of Eq. A. 22 are 
equal. Using the relation^ 4 , 









the desired inversion is obtained. With \) =1 and 

ARG = (*X» v/ t^-T^', Eq. A. 21 may be rewritten as 

ri n 

f(t') - ^ 

msi v / 2/^*1 



(A. 23) 



(A. 22+) 
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2 . 



- \ f c eKpL~ T ^ v ^'j 
^ s v/'ssw 



(A. 25) 



With V q = 1/s, Eq. A. 25 becomes 



p - ex^^T \|sVu)>] 



S^VsVoV 



(A. 26) 



The identity 



06 



4- - a E f. 



-hi 






L\fi^P 



+K 



(A. 27) 



is used to generate the desired inverse transform. Squaring 
Eq. A. 27 gives 

if = H Z M 

V>-o 




Letting x = 00 n /s in Eq. A. 28 and rearranging yields 

°° O(v-ti) 



(&) - 4 £ 



(A. 29) 



Solving Eq. A. 29 for l/s and substituting the value obtained 
into Eq. A. 26 results in 



\ZsV6J n ^’ 



(A. 30) 
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The inverse transform of Eq. A .30 is available 



( 8 ) 



There- 



fore with + 1 = tY) , 




(A. 31) 



A-3 : Inverse Laplace Transforms for a Gated Sine-Wave 

Input in Velocity in an Isovelocity Layer . 

The following acoustic quantities could not be 
completely solved. Although the expressions listed are 
incomplete solutions, it is worthy of note that the leading 
terms are the corresponding terms of the duct solutions 



the inverse Laplace transform of the quantity in brackets. 



in Section A-l. 




denotes 
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(A. 34) 




+■ 3c. 



u 



8 CL 






£Kf> [-T\/sW 






where ARG = 




7 



'Urr*- u(i-rY t Z's/A 





Loramel Function of two variables 



(14,15) 
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APPENDIX B 
EXPERIMENTAL DATA 



DISPLACEMENT 


displacement 


VELOCITY 


AXIS 


CROSSINGS 


AMPLITUDES 


AXIS 


CROSSINGS 


TPR I 


IN USEC 


(NORMALIZED) 


TPR I 


IN USEC 



STEP INPUT RUN S3 
C/AWN=0 .35 R=61.5CM 


A=2.54CM 


D=2 . 80CM TMP=20DEGC 


0.0 


0.175 


0.95 


1.10 


-0.379 


2.35 


2.85 


0.442 


4.35 


5.50 


-0.592 


7.65 


9.00 


0.660 


11.3 


13. 3 


-0.718 


16. 5 


18.5 


0.864 


21.7 


24.3 


-0.922 


28.2 


30.8 


1.000 


34.3 


38.0 


-0.898 


42.4 


45.9 


0.893 


50.3 


54.2 


-0.718 


59.2 


63.8 


0.718 


68.2 


73.5 


-0.592 


79.5 


84.0 


0.563 


91.0 


96.0 







STEP INPUT RUN S12 

C/AWN=0 • 52 R=61.5CM A=2.54CM D=4.15CM TMP=20DEGC 



0.0 

2.50 

5.10 

10.7 

17.7 
26 . 1 

36.6 

48.4 

61.3 

76.0 

91.7 

108.5 

126.0 

145.0 

164.0 

184.5 



0.198 

-0.396 

0.636 

-0.833 

1.000 

-0.972 

0;905 

-0.860 

0.838 

-0.838 

0.721 

-0.694 

0.622 

-0.613 

0.532 

-0.477 



NOT RECORDED 
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DISPLACEMENT 


displacement 


VELOCITY 


AXIS CROSSINGS 


AMPLITUDES 


AXIS CROSSINGS 


tpri in usec 


(NORMALIZED) 


TPRI IN USEC 


STEP INPUT RUN S6 






C/AWN=0.69 R = 6 1 


• 5CM A=1.Q0CM 


D = 5» 50CM TMP=20DEGC 


0.0 


0.414 


2.20 


3.50 


-0.687 


6.70 


9.60 


0.990 


14.4 


18.4 


-1.000 


23.2 


29.8 


0.919 


38*6 


45.8 


-0.909 


52.5 


64.0 


0.828 


73.5 


85.0 


-0.854 


93.0 


106.0 


0.798 


119.0 


131.0 


-0.727 


140.0 


168.0 


0.768 


169.0 


183.0 


-0.687 


193.0 


212.0 


0.687 




240.0 


-0.657 




268.0 


0.667 




298.0 


-0.581 




329.0 


0.616 




360.0 







STEP INPUT RUN S13 

C/AWN=0.85 R=61.5CM A=2.54CM D=6.78CM TMP=21DEGC 



0.0 
A. 20 

13.6 
26.3 

44. 6 
65.9 
90.2 
118.0 

147.5 

179.5 

212.5 



0.357 

-0.697 

0.871 

-0.874 

0.994 

-0.931 

0.777 

-0.834 

1.000 

-0.943 

0.783 



NOT RECORDED 
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DISPLACEMENT 
AXIS CROSSINGS 
TPR I IN USEC 



displacement 

AMPLITUDES 

(NORMALIZED) 



VELOCITY 
AXIS CROSSINGS 
TPR I IN USEC 



STEP INPUT RUN S1A 
C/AWN-1.00 R = 61 • 5CM A= 



0.0 


0 


7.00 


-0 


19.0 


0 


36.0 


-0 


59.5 


0 


89.2 


-0 


120.0 


1 


156.0 


-0 


190.0 


0 


235.0 


-0 



STEP INPUT 


RUN SA 


C/AWN-1 .00 


R=61.5CM A 


0.0 


0 


1.25 


-0 


3.A5 


0 


6. AO 


-0 


10.5 


0 


16.0 


-0 


21.9 


0 


28.2 


-1 


35.7 


1 


A5.2 


-0 


53.5 


0 


63.5 


-0 


7A.0 


0 


85.0 


-0 


98.5 


0 



D = 7 • 97CM TMP=21DEGC 

NOT RECORDED 



D-3.1ACM TMP-20DEGC 

0.80 
2.AA 
A. 90 
8.75 
13.3 
28.9 
25.2 
32. A 
AO. 5 
A9.2 

58.0 

68.5 

79.0 

91.5 



2.5ACM 

538 

813 

869 

875 

9A3 

956 

000 

863 

769 

619 

1 .OOCM 

176 

A12 

559 

598 

696 

81 A 

961 

000 

000 

971 

922 

863 

78A 

706 

667 
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DISPLACEMENT DISPLACEMENT 


VELOCITY 


AXIS CROSSINGS AMPLITUDES 


AXIS CROSSINGS 


TPRI IN USEC (NORMALIZED) 


TPRI IN USEC 


STEP INPUT 


RUN S9 






C/ AWN= 1.50 


R = 61 .5CM 


A= 1 ♦ OOCM 


D=4 • 7 1 CM TMP=2 1DEGC 


0.0 




0.307 


1.30 


2.30 




-0.586 


4.20 


6.50 




0.793 


9.20 


12.8 




-0.923 


16.5 


20.6 




1.000 


25.0 


30.6 




-0.964 


35.7 


42. 5 




0.979 


49.0 


57.0 




-0.943 


65.0 


71.5 




0.821 


79.0 


88.0 




-0.686 


98.0 


107.0 




0.750 


116.0 


125.0 




-0.700 


135.0 


145.0 




0.593 


154.0 


169.0 




-0.543 


178.0 


191.0 




0.521 


202.0 


217.0 




-0.543 


225.0 


239.0 




0.564 


246.0 


262.0 




-0.479 


269.0 


286.0 




0.486 


297.0 


312.0 








STEP INPUT 


RUN S5 






C/AWN=2.00 


R=61 • 5CM 


A= 1 • OOCM 


D=6.28CM TMP=20DEGC 



0.0 


0.351 


2.10 


3.80 


-0.649 


7.60 


11.3 


1.000 


15.5 


21.1 


-1.000 


26.2 


34.4 


0.979 


41.0 


52.5 


-0.918 


59.5 


72.5 


0.907 


79.5 


95.5 


-0.701 


101.0 


119.0 


0.731 


127.0 


145.0 


-0.639 


155.0 


176.0 


0.639 




203.0 


-0.639 




235.0 


0.639 




265.0 


-0.598 




295.0 


0.495 
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DISPLACEMENT 


displacement 


velocity 


AXIS 


CROSSINGS 


AMPLITUDES 


axis 


CROSSINGS 


TPRI 


IN USEC 


(NORMALIZED) 


TPRI 


IN USEC 



STEP INPUT RUN S7 
C/AWN=3.00 R=61.5CM 


A* 1 • OOCM 


D=9. 42CM TMP=20DEGC 


0.0 


0.432 


4.50 


7.50 


-0.892 


15.0 


21.0 


1.000 


31.0 


40.5 


-0.826 


53.0 


67.0 


0.822 


84.5 


100.0 


-0.784 


115.0 


132.0 


0.526 


150.0 


184.0 


-0.427 


210.0 


222.0 


0.376 




272.0 


-0.319 




312.0 


0.639 





88 



VELOCITY 
AXIS CROSSINGS 
TPR I IN USEC 



VELOCITY 

AMPLITUDES 

(NORMALIZED) 



GATED SINE INPUT RUN G1 W=583.3K.HZ 

C/AWN* 1 .00 R=6l • 5CM A=2.54CM D=7.97CM TMP=20DEG 



3.80 

9.60 

15.9 

21.5 
26.3 

31.9 

37.0 
42.2 

47.9 

53.0 

58.5 

63.5 

69.0 

75.0 

79.5 

84.5 

89.5 

95.0 

100.5 

105.0 

115.5 

120.5 

125.5 

131.0 

136.5 

141.5 

147.0 

152.5 

157.5 

164.0 

172.5 



0.279 

-0.669 

1.000 

-0.961 

0.747 

-0.870 

0.851 

-0.792 

0.909 

-0.870 

0.857 

-0.929 

0.864 

-0.909 

0.870 

-0.877 

0.922 

-0.864 

0.916 

-0.877 

0.890 

-0.883 

0.857 

-0.896 

0.864 

-0.877 

0.883 

-0.857 

0.896 

-0.831 
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APPENDIX C PROGRAM STEPINS2 
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N4=NUMBER OF EXPERIMENTAL MAX I MA ( POS I T I VE OR 
NEGATIVE) IN DISPLACEMENT. 



ampxz=experimentally determined values of peak 

AMPL I TUDES ( POS I T I VE OR NEGATIVE) IN 
PARTICLE DISPLACEMENT. 
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VNCOR=CORRECTION term IN VNRAT. SEE STATEMENT 601 
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1 DELT+TPRI ( 1 - 1 ) 

102 CONTINUE 
LO =22 

TEMP = TVAL ( I ) 

11 CALL BES(LO»ARG( I ) * 0 * BRESUL ( I ) ,S ) 



BSUM=TEMP*S(2) 
TEMP2=TEMP**2 
12 DO 13 L=4,L0»2 
TEMP=TEMP*TEMP2 
BSUM=BSUM+TEMP*S ( L ) 



r- 



















X 




Qd 




















z 


4C 0 


Qd 




O 




















M I 


~ ♦— 


UJ 






















4c 


Jt 


CM 


I- 




z 




















O 


\ r-H 






»— » 




















3 




O 




< 




















-I 


•-H M 


z 




X • 


















•— 1 


O 


1 X 


0 


• 


z 0 


















W- 


3 


X - 


U UJ 


isj UJ 


















1 — 




04 X 


UJ 


Qd 


Qd 


















< 


UJ 


— 3 


CO 0 


Z (X 


















X 


X 


w to 


#> 


3 


~ 3 


















Z 


I— 


4c 4c 


z 


UJ 


U 


















> 




4c — 




CD 


0 u 


















w 


X 


— Qd 


< 




UJ 0 


















u. 


0 


^ \ 


X 


to 


Qd 


















cO 


X 


1 - < 


z 


< 


Qd 1- 


















CO 


3 


CL 1 


isj 




3 — 


















< 




U H 




UJ 


U 


















+ 


Qd 


+ ~ O 


cO 


X 


U X 


















— * 


O 


I- 4* Qd 


t— < 


< 


O U 




















U 


w ^ < 




to 


*— 4 


















r-H 


z 


\ - w 


X 




to X 


















1 


ISj 


— z X 


Qd 


to 


< 5 


















•— » 




h 


UJ 


1 — I 


X 


















w. 


<3 


CL 4c 


I- 




1 — 


















1 “ 


z 


U < CD 




0 


0 < 


















< 


< 


1 — 3 


I- 


Qd 


z 


















Qd 




h \ l/) 


to < 


•— » *— » 


















Z 


z 


- U 3 


Qd 




to X 


















> 


>— 1 


— — 4c 


1— 1 


O 


to CL 


















w 


< 


4c 4c X 


Ll 


z 


O *- 


















LL 


X 


— 4- 4c 




< 


X 






to 


<— •* 










to 


z 


H S 4( 


UJ 


•> 


U to 








x 










CD 


ISj 


I CD ^ 


X 


Qd 


UJ 






— 


\ 










< 




X 1 ~ 


1 — 


O 


to z 








< 






• — • 






1— 


CM — I- 




U 


1— 1 »—l 






'-T 


1 






( — 1 




\ 


< 


— O CL 


UJ 


z 


X X 






3 


• 






1 






(X 


— Qd U 


Qd 


ISJ 


< X 






3 


H 




«*-> 


»— » 






z 


GO < + 


UJ 




UJ 






to 


w- 




X 






i“H 


ISj 


— 1— 


X 


to 


z *- 






3 


★ 




0 


1 — 


CM 


1 




O ~ — 


3C 




< UJ 






X 


z 


2 : 


U 


< 


4* 




CO 


►— *“H X 






0 






U 


3 


3 


z 


X 


H 




UJ 


1 — 






X 






» 


< 


to 


> 


z 


•» 


1 — 


1- 


^ X *- 






M h- 






O 


u 


CO 


1 


> 


CM < 


< 


II CM X 






t-H 








★ 


4c 


— * 


4c 


Qd 


3 


X ~ U 






to ^ 






*— 


• z 


l~» 




r-H 


z 


_l 


— CD | 






UJ 0 






»—» 


r-H 


> 


-w 


1— 1 




> 


< 


X 3 h- 






Z to 








★ 


U 


z 


W 


!— 1 




> 


3 CO — 






*— • 






O CM 


• 


< 


1— 1 


h- 


3 ~ 


UJ 


tO 3 ~ 






X X 






X — CM X 


4* U- 


< 


< 


r-H CO r-H 




II 






X •— 






< CNJ CM * X 


\ x 


x 


X 




CD 1 


z 








UJ 






* ~ 4c O LU 


• 


0 


z 


Z 


X < •— 


0 


Z 






h- • 




•— 1 


O to 4c 21 I— 


co u 


> 


> 


> 


II — 


— z 








x x 




w 


X 4c CL + 4c 


II 


II 


II 


II 


to 


■ — ■ •— » 


1“ 0 


s 






0 X 




3 


— X X <f X 


LU Z 




• — ■ 


Qd 


O X Qd UJ 


Qd ~ 


CM 






UJ 




< 


tO X 3 II X 


3 > 


»— * 


»— 1 


> 


Qd X X 3 


O CO 


w 






0 h- 




> 


3 3 1- X 3 


Z U 


w 


«*-* 


tO U 


— 1 - z 


CL CO 


X 






co 




I- 


CD K II 1 — 


< X 


1 - 


O 


X Qd + *—* 


UJ 


h- 






3 I- 


xO 


II 


II CM 0- II 


*- U- o < 


Qd 


<>*-♦- 


CO Qd 


< 






< < 


*— 1 X 


3 2E X *-h X 


z a u or 


U 




O -1 z 


~ Q_ 


Qd 






X 


II 


X 


3 3 X X 


o o z 


z 


X 


Ll_ 


u tu 0 


X X 


z 






h- z 


O UJ 


< CO 3 O 3 


uu > 


> 


< 


*— • 


< £> U 


f — UJ 


rsj 






►-* ISJ 


X 


I- 


U W H a h 



co 




co 



m 






UUUUUUUU<3<3UUU<3U 



97 



TEMP=TEMP*TEMP2 



€SUM=CSUM+(M-1 )*TEMP*S(M) 

17 CONTINUE 

ZNMA I N ( I ) =CSUM 

AXCROSZM=ZNMA I N ( I ) *ZNMAIN ( 1-1 ) 
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2130 IF(ZNRAT( I-1)-ZNRAT( I ) ) 2140,2160*2160 
2140 IF(ISIGN) 2180*2150*2190 
2150 IF(I-l) 2190 *2190 *2155 
2155 I S I GN=+ 1 

GO TO 2190 



2160 IF(ISIGN) 2190*2170.2180 
2170 IF(I-l) 2190.2190.2175 
2175 I S I GN=- 1 

GO TO 2190 
2180 MAXCT=MAXCT+1 
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PRINT 24.R.A,EN.D.C.TN*D0VERA,WN ♦ CAWN .CORFACPN ♦ CORF AC VN .CORFACZN 
24 FORMAT ( 4H1R =.E10.5»4H A »,El0.5*5H EN ».E10.5.4H D =.E10.5, 
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READ AO, IT 

CALL DRAW(NN,TPRI , VNRAT , 0 , 0 , 4h , I T , 0 , 0 , 0 * 0 , 0 , 0 , 9 ♦ 15 , 0 , LAST ) 

READ 40, IT 

CALL DRAW(NN»TPRI »VNMAIN»0,0»4H , I T ♦ 0 , 0 , 0 ♦ 0 ♦ 0 ♦ 0 , 9 , 1 5 , 0 , LAST ) 

READ 40, IT 
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tpri=delayed time. the time which has elapsed at a position r 

SINCE THE TRANSIENT FIRST REACHED THAT POINT. TPRI WILL RUN 
FROM 0 TO TPMAX IN INCREMENTS OF DELT. ‘ 
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DO 210 M=MSTART »N0 »2 
TEMP=-TEMP*TEMP2 
VAL=VAL+TEMP*B( M) 
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